arXiv:1502.01021vl [astro-ph.HE] 3 Feb 2015 


Mon. Not. R. Astron. Soc. 000.[T]-?? (2013) 


Printed 5 February 2015 (MN MgX style file v2.2) 


Relativistic jets shine through shocks or magnetic reconnection? 


Lorenzo Sironi/* Maria Petropoulou^t and Dimitrios Giannios^ 

^ Hm-vard-Smithsonian Center for Astrophysics, 60 Garden Street, MA 02138, USA 

^Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, IN 47907, USA 


Received / Accepted 


ABSTRACT 

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet 
flow is characterized by a high radiative efficiency and that the dissipative mechanism must 
be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have 
long been considered as possible candidates for powering the jet emission. Recent progress 
via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical 
grounds. We show that shock models are unlikely to account for the jet emission. In fact, 
when shocks are efficient at dissipating energy, they typically do not accelerate particles far 
beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection 
can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the en¬ 
ergy density of the magnetic field in the bulk flow is larger than the rest mass energy density. 
The emitting region, i.e., the reconnection downstream, is characterized by a rough energy 
equipartition between magnetic fields and radiating particles, which naturally accounts for a 
commonly observed property of blazar jets. 

Key words: acceleration of particles — galaxies: jets — gamma-ray burst: general — mag¬ 
netic reconnection — radiation mechanisms: non-thermal — shock waves 


1 INTRODUCTION 

Relativistic jets are ubiquitous in the Universe, with gamma-ray 
bursts (GRBs) and active galactic nuclei (AGNs) being two repre¬ 
sentative examples of high-energy sources powered by jets. Despite 
decades of research, the issue of what powers relativistic astrophys- 
ical jets is still unresolved. On the one hand, there is a strong theo¬ 
retical motivation t o consider them as magnetically-dominated ob¬ 
jects at their base ( iBlandford & Znai^ll977l : lBlandford & Pavnd 
Il982h . The strong magnetic fields threading a rotating compact ob¬ 
ject or the associated accretion disk serve to convert the rotational 
energy of the central engine into the power of the outflow. Part of 
the magnetic energy is used to accelerate the flow to relativistic 
speeds but, generally, magnetohydrodynamic (MHD) models pre¬ 
dict that the largest fraction of the energy remains locked in the 


magnetic field. 

.e., the jet arrives Poyntin 

g-dominated at the dissi- 

pation distance i 

Tchekhovskov et al.l2009 

ilKomissarov et alj2009l: 

lLvubarskvll200S 

). On the other hand, leptonic models of the radia- 


tive signature of blazars - a subclass of AGNs whose jets point 
along our line of sight - indicate that the energy densities of mag¬ 
netic field (7 b and radiating particles (7e at the e mission region do 
not differ by more than one order of magnitude jReadheadll 19941 : 
ICelotti & Ghisellinil2008l:lGbisellini et'^l2014t) . 

Our ignorance on the source of jet power allows for 
two very different scenarios regarding the dissipative mecha- 
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nism behind the jet emission. Hydrodynamical flows can dissi¬ 
pate energy and efficiently accelerate particles at s hock fronts 
(e.g.. lHeavens & Drury||l988l : lRees & Meszaroslll994l) . whereas in 
strongl y magnetized flows the ef ficiency of shocks is greatly re¬ 
duced jKennel & Coronitj|l984hh . In this case, a more likely can¬ 
didate for p owering the jet em i ssion is the process of mag netic re¬ 
connection dSpruit et al.ll200ll : IPrenkhahn & Sprui3l200A . where 
the annihilation of field lines of opposite polarity transfers the field 
energy to the particles. 

Observations of relativistic jets can be used to infer the prop¬ 
erties of the emitting region. Yet, these do not directly represent the 
typical conditions in the bulk flow. Emitting regions are, by defi¬ 
nition, special places were intense energy dissipation has occurred. 
Both the shock downstream and the reconnection outflow qualify 
as such regions. On the contrary, the upstream medium (in the case 
of shock dissipation) or the inflow region (for the case of recon¬ 
nection) describe more closely the bulk of the jet. In this sense, the 
observational inference on the ratio (7 b/( 7c, applies directly to the 
emitting region and not to the bulk of the jet flow. Hereafter, we 
identify the (shock or reconnection) upstream with the large-scale 
jet and the downstream with the emitting region. 

The properties of the emitting region are related to those of the 
bulk of the jet in a way that depends on the dissipative process. In 
shocks the magnetic field is generally compressed and its strength 
is, thus, increased: the magnetic energy per particle is larger in the 
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downstream than in the upstreamQ In contrast to shocks, magnetic 
reconnection dissipates magnetic energy: the downstream (i.e., the 
outflow) is less magnetized than the upstream (i.e., the inflow). In 
fact, the more efficient is reconnection in accelerating particles, the 
more particle-dominated the downstream region becomes. 

In both shocks and reconnection, the partition of energy 
between magnetic fields and particles can be described self- 
consistently via fully-kinetic particle-in-cell (PIC) simulations. Re¬ 
cent PIC progress has been remarkable. First-order Fermi accel¬ 
eration in relativistic shocks has been demonstrated from first 


princ iples i Snitkovsk^ 2008al l5: iMartins et all l2009l : iHaugbdll^ 
1201 ll: ISironi etalJ 1201^ However, as expected from theoreti¬ 
cal arguments jBegelman & Kirlj|l990ll . no particle acceleration 
takes place when relativistic shocks are quasi-perpendicular, i.e., 
where the field lin es are nearly orthogonal to the shock direc- 
tion of propagation jGallant et alJl993 : ISironi & SDitkovskvl2009L 
l201lh . Recent PIC simulations have also demonstrated that mag¬ 
netic reconnection in magnetically-dominated plasmas is a fast, 
unsteady process that leads to efficient parti cle acceleration (e.g., 
ISironi & SDitkovskvl2014l : lKagan et alj2015l . for a review). In both 
shocks and reconnection, leptons pick up a large fraction of the dis¬ 
sipated energy, so electrons are the most likely candidates for pow¬ 
ering the observed emission (i.e., supporting the so-called leptonic 
scenario). 

In this work, we use jet observations to probe the nature of the 
dissipative mechanism and the properties of the bulk of the jet. For 
this, we exploit the latest PIC findings about shocks and magnetic 
reconnection as well as three basic facts about jet emission: (i) jets 
are efficient emitters, (ii) the radiating particles have extended, non- 
thermal distributions, and (iii) the energy densities stored in the 
magnetic field and in radiating particles do not differ more than one 
order of magnitude^ The paper is structured as follows. In Section 
2 we summarize these basic observational constraints. In Section 3 
we apply them to two alternative dissipative scenarios: shocks and 
magnetic reconnection. We conclude in Section 4. 


2 INFERENCES EROM JET OBSERVATIONS 

There is a consensus that relativistic jets are efficient emit- 


power into gamma-ra 

y (prompt) emission (e.g.. Berger et alj200C 

Panaitescu & Kuman 

2 OO 2 I: lYost et all I 2 OO 3 I: ICranot et aljl2006 

Nava et al.l 2013h. Blazars have similarly high radiative efficien- 

cies. Recently. Ghisellini et al.|(|20l4l modelled the spectral energy 


distribution (SED) of a large sample of blazars and, assuming one 
proton per radiati ng electron, they inferred a radiatiye efficiency 
/ ~ 0.1 (see also lCelotti & Ghisellinill^OSh . The efficiency may 
be even higher, if blazar jet s are pair-loaded, nam ely they contain 
more leptons than protons l lGhisellini et al.ll2014h . It follows that 
the dissipative mechanism behind the jet emission must be able to 
transfer a large amount of energy to the radiating particles, achiev¬ 
ing a high radiative efficiency (/ > 0.1). Given that / is already 


high, the dissipation process must involve a large fraction of the jet 
fluid, in order to explain the observed photon luminosities without 
pushing the jet energetics to extreme values. 

Multi-wavelength jet observations typically suggest that the 
radiating particles are accelerated into a non-thermal power-law 
distribution. The blazar SED has a characteristic double-humped 
shape (e.g.. lUlrich et al.iri997l : lEossati et alJll998l) . with a broad 
low-energy component extending from the radio up to the UV band, 
and in some extreme cases up to > 1 keV X-rays l lCostamante et al.l 
l200lh . The high-energy component extends across the X-ray and 
y-ray bands, with peak ene rgy around 0.1 T eV, although this is 
not always clear (see, e.g., lAbdo et ^l201ll for Mrk 421). The 
low-energy hump is believed to result from synchrotron emis¬ 
sion from relativistic electrons. Its spectral shape requires, in gen- 
eral, a power-law (or bro ken power-law) particle spectrum (e.g., 
ICelotti & Ghisellinill2008ll . which implies non-thermal particle ac¬ 
celeration. Furthermore, the high-energy hump of blazars, which in 
leptonic scenarios is explained as inverse Compton emis sion, ex¬ 
tends in many sources up to 100 GeV or higher (e.g., see lAlek^ 
I 20 II for PG 1553-1-113). This also points to electron acceleration 
up to Lorentz factors ye <: 10^, or even higher for TeV-emitting 
blazars The same holds true in synchr otron models for the prompt 
emission of GRBs jPaigne et alj|201 ll . and references therein)01n 
summary, the dissipative mechanism that operates in jets has to be 
able to accelerate particles well beyond their thermal energies. 

In blazars both the synchrotron and the Compton humps are 


visible and their relative brightness constrains the ratio of energ y 
densities UbIUc at the emitting region (e.g.. ISikora et alj|2009l) . 
The fact that in BE Lacs the two bumps have comparable luminos¬ 
ity indicates that, in th e synchrotron self-Compton (SSC) model, 
UbIUc ~ 1 (see also iMastichiadis & Kirk|[l9^ . However, for 
gamma-ray luminous blazars, such as flat spectrum radio quasars 
(FSRQs), where the high-energy component is much more lumi¬ 
nous than the low-energy one, the SSC scenario would predict that 
UbIUc <k 1, i.e., far from equipartition. Yet, a rough equiparti- 
tion between the magnetic field and radiating electrons is still in¬ 
ferred in such blazars, in the ext ernal Compton (E C) interpreta¬ 
tion of their gamma-ray emission l lSikora et al.l2009l) . In summary, 
detailed modeling of many blazars within the SSC or EC scenar¬ 
ios points to the fact that, quite generally, UbIUc ~ 1, but some 
emitting regi ons may be moderately magnetically-dominated, i.e., 
Ub/Uc ~ 3 dChisellini et alj|2014l) . The rough equipartition be¬ 
tween emitting particles and magnetic field in the blazar emitting 
region is also suggested by th e observed surface brightness tem¬ 
peratures at radio wavelengths 1 R£adheacijV94|: lL^teenm^ et al.l 

.1120131). 


119991 : iHoman et al.ll2006l : iHovatta et al. 120131) . In this work, we 
adopt the range 0.3 < f/B/l/e < 3 as represent ative of the physical 
conditions in the emission region of blazars 1 Cop pj_& Aharonia^ 


ll999l:ICelotti & Chisellinil2008l:IChisellini et alj2014h. In the ( 

of GRBs, no high-energy bump in the SED can be securely asso¬ 
ciated with the prompt emission, making the ratio at the 

emitting region harder to quantify. 

In summary, any model for the jet emission has to account for 
the fact that both blazars and GRBs are characterized by (i) a large 


^ This condition is verified in quasi-perpendiculai* shocks (i.e., where the 
field lines ai‘e nearly orthogonal to the shock direction of propagation), 
which is the most relevant field configuration for relativistic jets, due to 
Lorentz transformation effects and shock compression. 

^ We remark that the observational inferences (i) and (iii) ai‘e model- 
dependent, relying on the assumption of a leptonic origin for the observed 
emission. In this study, our working hypothesis is that the radiating particles 
are leptons, and our results should be interpreted in this framework. 


^ In hadronic models for the gamma-ray emission from blazars, the re- 
quired maximum Lorentz factor of protons may be even mo re extreme 
iAhai'oni^l2QQ0l : iBottcher et al.l2of3;IPetropoulou et alj2Q15h . 

^ Note however that the GRB prompt signa ture can be the result of 
photospheric emission from ther mal electrons iGhisellini & 0010151119991 : 
lGianniosl2QQdlPe’er et al.ll2QQ6h . 
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radiative efficiency / > 0.1 and by (ii) extended non-thermal distri¬ 
butions of the radiating particles. In addition, blazars are also char¬ 
acterized by (iii) a rough equipartition of energy density between 
radiating electrons and magnetic field. 


3 CONSTRAINTS ON THE DISSIPATIVE MECHANISM 


In this section we investigate whether the observational inferences 
- on the dissipative efficiency, the equipartition between particle 
and magnetic energies, and the requirement of non-thermal particle 
acceleration - can constrain the dissipation mechanism responsible 
for the jet emission. Shocks and magnetic reconnection are the most 
commonly invoked dissipation processes. If the jet is dominated by 
kinetic energy flux at the dissipation distance (low-cr flow)0 then 
shocks are natural candidates for the required dissipation. On the 
other hand, if the jet remains Poynting flux dominated at the dissi¬ 
pation distance, magnetic reconnection is a more plausible candi¬ 
date for powering the je t emission, since shock dissipation is sup¬ 
press e d in high-(T flows ^Kennel & Coronitiiri984bl : lOiannios et al.l 
l2008l: iNaravan et alj|201 ih . In this case, however, a process that 
triggers the reconnection of magnetic field lines is required, e.g., 
MHD instabilities at the dissipation radius or large-scale changes 
in the magnetic field polarity of the flow i Romanova^^ovela^ 


1V9^ Spmit et ffi]200l|: Drenkhahn & Spruit|2002: LYUtiko^^200f 


Gianniosll2010l;IZhans & YaiJl201ll;lMcKinnev & Uzdenskvi2012; 

IparfrevLljS. 


To minimize the energetic requirements of the dissipation 
mechanism, we assume, for both shocks and reconnection, that rel¬ 
ativistic electrons » 1) are responsible for the observed ra¬ 
diation and that they lose most of their energy radiatively. In what 
follows, therefore, the terms radiative and dissipative efficiency will 
be used interchangeably. We also assume that electrons pick up a 
large fraction of the dissipated energy in both models, while the 
rest of the dissipated energy goes into protons, or ions more gener¬ 
ally. Thus, the electron energy density Uc in the emitting region is 
parametrized as a fraction of the overall energy density e, i.e., Uc = 
e, where e does not include the rest mass energy density. For the 
discussion that follows, we set = 0.5 for electron-proton flows 
and = 1 for electron-positron flows. These values are physi¬ 
cally m otivated on the basis of PIC simulations, for both relativistic 
shock s l lSDitkovskvll2008j : ISironi & Spitkovskvll201 lUSironi et al.l 


2013h and 


magnetic reconnection i 


Melzani et al.ll2014l) . 


(e.g., Sironi & Spitkovskvll2014l: 


3.1 Shocks 

In the case of shocks, the downstream (emitting region) and 
upstream (large-scale jet) properties can be directly connected 
through the jump conditions at the shock. Let us first consider a 
hydrodynamic shock that corresponds to a non-magnetized jet. For 
a cold upstream medium, the kinetic energy per unit rest mass in 
the downstream region is found to be e/pc^ = yrei “ 1. where 

7rci = 1/ ^1 - and j3„;i = Vrei/c is the relative velocity between 
the upstream and downstream regions, in units of the speed of light. 


^ Throughout this work, the magnetization of the fluid is defined as cr = 
where both the magnetic field strength B and the mass density p 
are measured in the rest frame of the fluid. 


The dissipative efficiency, which, under our assumptions, is equal 
to the radiative one, is then 


u = 


e + pc^ 



( 1 ) 


which in the ultra-relativistic regime (yrei 1) reduces to /sh ~ fe- 
For our physically-motivated choice of = 0.5 - 1, the above 
relation yields an efficiency in the range of the observed values. 

In general, the upstream region is magnetized, and in this 
case the efficiency drops, unless 7^1 increases accordingly. Thus, 
shocks formed in a flow with any appreciable magnetization must 
be mildly relativistic or relativist ic, in order to satisfy the efficiency 
requirement /si, >0.1 (see, e.g., iMimica et al.ll2009l) . This will be 
shown with detailed numerical examples at the end of this section. 
Similar to the hydrodynamic case, the downstream and upstream 
quantities can be related through the MHD jump conditions (for a 
treat ment of the MHD jump conditions in relativistic sho cks see, 
e.g.. iKennel & CoronitJ[l984al : IAddI & CamenzindfTgSSl) . In the 
following, we introduce the notation g,, with i = 1 for upstream and 
i = 2 for downstream, to refer to quantities of region i measured in 
the respective rest frame. The quantities to be determined and tested 
against the observational constraints are the efficiency (/i,) and the 
ratio of the magnetic to particle energy densities (UBalU^a) ss a 
function of the upstream magnetization cr|. From this point on, we 
will refer to cti simply as cr, in order to use a uniform notation 
throughout the text. The efficiency is defined as 


/sh - 


U, 


e,2 


62 + P2C~ + Db,2 


( 2 ) 


where (/e ,2 = /e 62 . Ub ,2 = By^n and p 2 is the proper density of the 
downstream region. The ratio of magnetic to electron energy is 


Ub,2 

~2 


F- 1 1 

/3p,2 


(3) 


where £2 = F 2 /(f - 1), f is the adiabatic index ranging from 5/3 
(non-relativistic ideal gas) to 4/3 (ultra-relativistic gas), and / 3 p ,2 = 
^ 2 /f^B ,2 is the plasma beta parameter for the downstream region, 
while P 2 is the gas pressure. 

The parameter /i, as defined above represents the fraction 
of post-shock energy that is deposited into electrons (or gener¬ 
ally, pairs), and thus available to be radiated. Yet, the observa¬ 
tional requirement of efficient dissipation (/i, > 0 . 1 ) is not the 
only constraint imposed on the candidate dissipation mechanism. 
Unless shocks can accelerate particles to energies well beyond the 
peak of their thermal distribution, dissipation at shocks is not a 
promising process for explaining the observed non-thermal emis¬ 
sion from jets (see also Section H. In fact, it is well known 
that Fermi accele ration cannot operate in the so-called superlumi¬ 
nal shocks (e.g.. [K irk & Heavend 1198^: iBege lman & Kl3ll99d : 


ISironi & Snitkovs^ 20091 201in^oni et al]2013 i. There, parti¬ 

cles moving along the magnetic field at the speed of light can¬ 
not outrun the shock, whose speed in the pre-shock medium is 
vi ~ c, for ultra-relativistic flows. Since the Fermi process requires 
repeated crossings of the shock, this implies that Fermi accelera¬ 
tion in superluminal relativistic shocks is extremely inefficient. For 
this reason, we also require the shock to be subluminal, or equiv- 
al ently the condition [i\ / cosSj < 1 should hold (see, e.g., eq. ( 1 ) 
in iKirk & Heavenslll989h . where Pi = Vi /c and 9i is the angle be¬ 
tween the shock normal and the magnetic field as measured in the 
upstream frame. 

We solve the MHD jump conditions for a cold upstream 
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mediurr@ and two representative shock obliquities, i.e. = 30° 
and 60°. The results are presented in Fig. [T] and our main conclu¬ 
sions can be summarized as follows: 


• For cr < 0.1, the efficiency does not strongly depend on the 
flow magnetization. 

• The shock must be mildly or ultra-relativistic in order to be 
efficient, since /^h > 0.1 only for yrdiSmi <: T which is consistent 
with the argument bv iRees & MeszarM dlQQdh . 

• The ratio of magnetic to particle energy densities in the down¬ 
stream increases with the upstream magnetization. For cr =: 0.001 - 
0 .1, there is a small range in yrd/Jrei (dependent on cr) leading to 
0.3 < UB, 2 lUe ,2 < 3. This conclusion does not depend significantly 
on the magnetic field inclination relative to the shock front. 

• For oblique shocks with angles 6i ;£ 30°, there is a small 
parameter space where all of the three constraints are satisfied, 
namely 0.03 < fsh < 0.3, 0.3 < Ub, 2 IUc ,2 ^ 3 and the shock is sub¬ 
luminal. This parameter regime involves mildly relativistic shocks 
withy6jei7rd - 1 and moderate magnetization cr ~ 0.1. 

• For substantial magnetic field inclination (di > 45°), there is 
a strong tension between the efficiency constraint, that favors fast 
shocks, and the subluminality requirement. For 0; = 60°, we find 
that the efficiency for subluminal shocks is at most ~ 3%. 


The analysis above assumes that the magnetic field in the 
downstream region is merely the result of shock compression (as 
opposed to fields generated via plasma instabilities at the shock 
front) and that the upstream region is sufficiently magnetized for the 
superluminal constraint to hold, i.e., cr > 10“^ i n electron-positro n 
shocks and cr > lO^"* in electron-proton shocks dSironi et al.l2013l) . 
We showed that the observational requirements of 0.03 < pt $ 0.3, 
0.3 < Usa/Uca i 3 and the subluminality constraint for efficient 
particle acceleration are simultaneously satisfied only for certain 
shock obliquities and for cr > 0.01. In this case, the downstream 
field is likely to be dominated by the shock-compressed compo¬ 
nent, in agreement with our assumptions. 

These assumptions, however, may not hold in other environ¬ 
ments, such as the external shocks of GRBs. There, the shock 
propagates into the interstellar medium, whose magnetization is 
extremely small (cr ~ 10“*). Modelling of the GRB afterglow 
emission implies t hat the emitting region is far from equipartition, 
havin g I/e » Ub i Kurnar & Bamiol Duraiill200^ . l20ld : iLemoinel 


l2013l:fsironi et alj|2013h . Thus. in the context of GRBs, our argU' 

ments are applicable only to the “internal” jet emission mechanism 
and not to the afterglow phase. 


3.2 Magnetic reconnection 


In reconnection, magnetic field energy is transformed into parti¬ 
cle energy. Assuming that half of the magnetic energy is converted 
into internal energy (quite a realistic assumption, as we demon¬ 
strate below with PIC simulations), and assuming fast cooling elec¬ 
trons, the expected radiative efficiency is /rec ~ 0.5 fe cr/{(T+2). For 
~ 0.5-1, as found in PIC simulations of electron-positron (yield¬ 
ing 1) and electron-proton (giving 0-fi) reconnection 


(e.g., Sironi & Sr)itkovsk^l2014l : lMelzani et al.ll20l4) . this implies 
that cr > 1 for reconnection to be energetically viable in powering 
the jet emission. Hence, we deal with relativistic reconnection, in 


^ We numerically solve eqs. (25), (28), (29), (31), (33) and (37) in 
IApdI & Camenzin3 jl988h . while we use eqs. (13)-(14) in lServic3 jl986l) 
for the dependence of the adiabatic index on the downstream temperature. 
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e,=6o“ 


-2.5 


-0.5 
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Figure 1. Plot of yrei/^rei as a function of the upstream magnetization cr 
for two different inclinations of the upstream magnetic field relative to the 
shock normal (measured in the upstream frame): = 30° (top panel) and 

60° (bottom panel). The black, red and green colored regions lead, respec¬ 
tively, to efficiencies 0.03 ^ Ah ^ 0.3, ratios of magnetic to electron energy 
densities 0.3 < I7B,2/I^e,2 3 and subluminal configurations. For param¬ 

eters lying in the purple-hatched region no fast magnetosonic shock solu¬ 
tions can be found. In both panels we assume = 0.5, as appropr iate for 
electron-proton flows jSironi & Spitkovsk4l201ll : ISironi et alj2013h . 


which the mean magnetic energy per particle is larger than the rest 
mass energy, or equivalently cr > 1. 

Recent progress in the study of relati vistic reconnection via 
PIC simulations has been remarkable (see iKagan et al.l [20151 . for 
a review), demonstrating that reconnection in high-cr plasma is 
fast and unsteady. The reconnection downstream is composed of 
plasmoids that contain the reconnected plasma. In two dimen¬ 
sions (2D), they appear as overdense magnetic islands (see the 
2D density pattern in the top panel of Fig. connected by thin 
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Figure 2. S tructure of the reconnection layer at tUpf = 3000, from a 2D simulation of cr = 10 anti-parallel reconnection in electron-positron plasmas, as 
described in ISironi & Soitkovskvl ilOldll . The box extends along x over ~ 6550 c/tUp (65536 cells), and along y over ~ 6000 c/tUp (~ 60000 cells), but we only 
show a subset of the domain, to emphasize the small-scale structures in the reconnection layer. Here, c/tUp = lAnn^e^ is the plasma skin depth, where 

fiQ is the particle number density far from the reconnection layer. We present (a) the particle number density in units of no, with overplotted magnetic field 
lines; (b) the particle density, with artificially colored in black the regions that do not meet the selection criterion y > 1.1 for the reconnection downstream; (c) 
the magnetic energy density normalized to the rest mass energy density far from the reconnection layer, eg = B^/8np(?, where p = mno is the mass density 
far from the reconnection layer; and (d) the kinetic energy density normalized to the rest mass energy density far from the reconnection layer. 


X-lines. The plasmoid instability further fragments each X-line 
into a series of smaller islands, separated by X-points (e.g., see 
the small islands at 700c/tUp < x < 1400c/tUp in the top 
panel of Fig. where c/tUp is the electron skin depth). At the 
X-points, the particles are not tied to the magnetic field lines 
and they get accelerated along the reconnection electric field. 
The late-time particle spectrum integrated over the whole recon¬ 
nection region is a power-law, whose slope is -2 for cr = 10 
(as commonly assumed in blazar SE P modeling), and becomes 


harder with increasing magnetization (Sironi & Soitkovskvl 

2014; 

see also. e.g.. IZenitani & Hoshind 12001; Jaroschek et al. 

2004; 

Bessho & Bhattachariee 

2 OI 2 I; Kaean et al. 2013: Cerutti et al. 

2014 GuoetalJ 2014; 

Werner etalj 2014). Efficient particle ac- 


celeration to non-thermal energies is a generic by-product of the 
long-term evolution of relativistic reconnection in both two and 
three dimensions. In three dimensions (3D), the so-calle d drift-kink 
mode corrugates the reconnectio n layer at early times toaughtonl 
1 19981 ; IZenitani & Hoshin31200 Sh . but the long-term evolution is 
controlled by the plasmoid instability, that facilitates efficient 
particle acceleration, in an alogy to the two-dimensional physics 
(ISironi & Spitkovskvll20l4) . 

In summary, recent PIC simulations of relativistic reconnec¬ 
tion have convincingly demonstrated that for cr > 1 (of interest for 
relativistic jets), the accelerated particles populate extended non- 
thermal power-laws. We now investigate whether relativistic re¬ 


connection can provide the required efficiency /kc and kinetic-to- 
magnetic energy ratio U^IUb inferred from jet observations. 

The dissipated magnetic energy in the reconnection down¬ 
stream is distributed between particles and magnetic fields. The 
plasmoids appear as inhomogeneous structures, with the core dom¬ 
inated by particle energy (panel (d) in Fig. O, whereas the outskirts 
are magnetically dominated (panel (c) in Fig. O. Because of the 
inhomogeneity of the plasmoids, it is meaningful to define the dis¬ 
sipative efficiency as a volume- or surface-averaged quantity in 3D 
and 2D, respectively, summed over many plasmoids: 

2,X,t/edF,- 

" ZjyXe+pc^ + Un)dV, 

where V, is the volume (or surface) of each plasmoid, 11^ is the 
kinetic energy of the emitting leptons (electrons alone for electron- 
proton reconnection, or both species for electron-positron recon¬ 
nection), while the sum runs over all the plasmoids in the recon¬ 
nection layer (see the colored regions in Fig. Ob)). 

The fact that the reconnection upstream is taken to be cold 
facilitates the identification of the regions where dissipation of en¬ 
ergy has taken place, as we now explain. We select the plasmoid 
volume Vi where to compute the integrals above (or actually sur¬ 
face, for 2D simulations as in Fig. O by selecting all the regions 
where the mean particle Lorentz factor is y > 1.1 (more precisely, 
in electron-ion reconnection we employ here the mean ion Lorentz 
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Figure 3. As a function of the flow magnetization, we plot the dissipa¬ 
tive efficiency (top panel) and the volume-averaged {UeHUe + Ub)) (bot¬ 
tom panel) for relativistic reconnection i n electron-positron plasmas, from 
a series of 2D simulations as described in isironi & Soitkovskvl J2Q14h . The 
black lines refer to the case of anti-parallel fields (so, BgjBo = 0, where 
Bg is the guide field and Bq is the alternating field), whereas the effect of 
the guide field is shown by the filled colored circles, as indicated in the leg¬ 
end. We note that the range of values 0.3 ^ UqIUb ^ 3 infeiTed from the 
observations con'esponds to 0.2 < f/e/((/e + Uq) < 0.75 for the case of 
homogeneous plasmoids. 


factor). The inflow region is cold and move s at the reconnection 
speed , which for cr > 10 is Vj-ec 0.15 c (e.g. JSironi & Spitkovskvl 
[ 2 on . so that the mean Lorentz factor there is only ~ 1.01, 
which is below the threshold we employ. In Fig.j^b), we only dis¬ 
play the regions where the condition 7 > 1.1 is met (otherwise, the 
color is artificially set as black), showing that this criterion provides 
an excellent identification of the reconnection downstream (i.e., of 
the magnetic islands). We have employed this criterion across the 
whole range of magnetizations we have investigated (cr = 1 - 50) 
and for different strengths of the guide field 6 g orthogonal to the 
annihilating fields Bq (from B^IBq = 0 up to 3), finding always an 
excellent spatial correlation with the reconnection downstream. In 
the following, we define cr = B\l‘XnpP, i.e., excluding the contri¬ 
bution of the guide field. 

We now explain how to extract the equipartition parameter 
UcIUb- Fits to observations are performed under the assumption 
of a homogeneous emitting region (the so-called one-zone homo¬ 
geneous emission models). They infer 0.3 < UcIUb :£ 3, which 
can also be expressed as 0.2 < U^I{U^ + Ug) < 0.75. Since the 
reconnection downstream is inhomogeneous (see the magnetic is¬ 
lands in Fig. 13, we need to define a measure of Uel{Uc + Ub) in the 
downstream region that closely corresponds to what the homoge¬ 
neous models probe. In the limit of fast cooling particles, we should 
calculate the volume-average (surface-average in 2D simulations) 


Table 1. Results from 2D PIC simulations in electron-positron plasmas 
without guide field, as a function of the flow magnetization. 


(T 

1 

3 

10 

20 

30 

40 

50 

free 

0.26 

0.44 

0.52 

0.55 

0.54 

0.53 

0.52 

{UJiU, + Ub)) 

0.75 

0.74 

0.72 

0.72 

0.71 

0.71 

0.70 


Table 2. Results from 2D PIC simulations in electron-positron plasmas, as 
a function of the flow magnetization and the guide field strength. 



BgIBo 

0 

0.1 

0.3 

0.6 

1 

3 

(T = i 

free 

0.26 

0.25 

0.21 

0.15 

0.11 

0.03 

(T = i 

{UJiUc + Ub)) 

0.75 

0.70 

0.57 

0.36 

0.22 

0.04 

0-= 10 

free 

0.52 

0.48 

0.36 

0.24 

0.13 

0.03 

0-= 10 

{UJ{U, + Ub)) 

0.72 

0.66 

0.52 

0.35 

0.22 

0.04 

o- = 50 

free 

0.52 

0.50 

0.35 

0.25 

0.14 

0.04 

o- = 50 

{UU(Ut + Fb)) 

0.70 

0.63 

0.49 

0.35 

0.21 

0.05 


of UcliUc + Ub) weighted with the electron energy density Ut, and 
then sum over all the plasmoids: 

/ De \ _ fvj Pc+Pb 

\U, + UbI~ Zjy^.U^dV, 

The reason for our choice of the weighting factor Ue is that re¬ 
gions with low electron energy density cannot contribute much to 
the emission. Although in regions where Ub is larger than U^ (e.g. 
in the outskirts of the plasmoids, compare Fig. [2lc) and (d)), the 
synchrotron cooling rate is significantly higher than in regions with 
lower Ub (e.g., the plasmoid cores), the total radiated power in the 
fast cooling regime is always limited by the one injected into the 
emitting particles. As a side note, we remark that the choice to in¬ 
tegrate UcKUc + Ub) rather than U^Ub in Eq. (|5) is motivated by 
the fact that the integral would otherwise be largely dominated by 
the very center of magnetic islands, where Ub 'K U^ (see Fig.|3c) 
and (d)). Our choice is then appropriate for comparing to observa¬ 
tional inferences based on homogeneous emission models. 

By integrating over the surface of the magnetic islands, we 
compute /rec and (U^KU^ + Ub)) as described above. Our main re¬ 
sults are shown in Fig.[3for the case of electron-positron reconnec¬ 
tion and in Fig.|4]for electron-proton reconnection. Our conclusions 
can be summarized as follows (see also Tables [T13 : 

• Jet observations do not impose strong constrains on the recon¬ 
nection model as long as the reconnection is relativistic. For cr > 1 
and guide fields weaker that the reconnecting field (Bg < Bq), all 
constraints are satisfied. We have checked that this conclusion holds 
for a broad range of conditions: pair and electron-ion plasma, and 
upstream magnetizations as large as cr ~ 50. 

• In electron-positron plasmas (top panel in Fig. |3, the recon¬ 
nection efficiency asymptotically approaches ~ 0.5 for cr > 10 
and no guide field. In electron-ion plasmas (colored lines in the top 
panel of Fig. |4), it is roughly half of that value, since in the nu¬ 
merator of Eq. ([ 4 ) only electrons contribute (as opposed to both 
species, for electron-positron plasmas). Eor lower cr, the rest mass 
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energy appreciably contributes to the denominator of Eq. which 
results in a decrease of the reconnection efficiency. 

• As a result of the efficient field dissipation, the energy density 
of the radiating particles is larger than that of the magnetic field at 
the emitting region by a factor of a few, both in electron-positron 
plasmas (black line in the bottom panel of Fig.^ and electron-ion 
plasmas (colored lines in the bottom panel of Fig.|4j. Once again, 
the difference between the colored lines and the black line in Fig.|4] 
is due to the fact that, in electron-positron reconnection (black line), 
both species contribute to the kinetic energy (7e of emitting parti¬ 
cles that enters in Eq. 

• For weak guide fields ( 6 g/So < 0.5), the magnetic energy 
is very efficiently dissipated (/jec <: 0.2 in electron-positron plas¬ 
mas). For progressively stronger guide fields, the radiative effi¬ 
ciency drops (see the colored circles in the top panel of Fig. [3} 
and the emitting regions turn magnetically dominated (colored cir¬ 
cles in the bottom panel). The transition happens gradually around 
Sg/Sfl ~ 0.5, independently of tr. 

• For electron-positron reconnection, we find that the ratio be¬ 

tween particle kinetic energy and magnetic energy has an upper 
limit of ~ 3, see Fig.[3lb) (and ~ 2.5 for electron- 

proton plasmas, see Fig. |4tb)). This is independent of a and is 
reached in the limit of weak guide fields. When increasing the 
guide field strength, the emission region gets more magnetized 
and can become arbitrarily small. Flowever, we remark that 

when the emission region is strongly magnetized, the dissipation 
process is inefficient (i.e., 1 ), so the resulting observational 

signature is very weak. It follows that in reconnection-dominated 
systems there is an observational bias to infer equipartition. 

• We can provide a convenient parameterization of the results 
presented in Fig. [3] and Fig. |4] (see also Tables [Tl^ . assuming that 
a fraction ~ 55% of the magnetic energy in alternating fields (so, 
excluding the guide field) is converted into particle kinetic energy, 
and that the guide field energy is not dissipated. We obtain 


/,ec 


( 6 ) 


where fe ~ 1 in electron-positron plasmas and ~ 0.5 in electron- 
proton plasmas. Similarly, the expected equipartition parameter is 


-^— = 1 3- ^—— - (7) 

(/e -I- (/b 0.55 fe + 0.45 -r (BJBof 

where the extra factor of 1.3 accounts for the fact that the integrand 
in Eq. has been weighted with the electron energy density t/c. 
These two formulae provide a satisfactory description of the trends 
observed in Fig.^and Fig.|4] 

• Our results for electron-ion reconnection are nearly indepen¬ 
dent of the numerical choice of the mass ratio, which is constrained 
in PIC simulations to be smaller than the realistic value, for compu¬ 
tational convenience. In fact, the three colored lines in Fig. |4] (blue 
for ;?7p/me = 6.25, green for m^jm^ = 25 and red for m^jm^ = 100) 
nearly overlap. It follows that our results can be confidently applied 
to the case of a realistic mass ratio. 


4 DISCUSSION AND CONCLUSIONS 

For a long time, theoretical modeling of relativistic astrophysical 
jets has faced a conundrum: if jets are launched as magnetically- 
dominated flows, as most theories predict, why do they appear to 
be in rough equipartition between magnetic fields and radiating par¬ 
ticles at the emission region? Does the Poynting flux convert into 



Figure 4. As a function of the flow magnetization, we plot the dissipative 
efficiency (top panel) and the volume-averaged {Uel{Ue + Ub)) (bottom 
panel) for relativistic reconnec tion in electron-ion plasmas , from a series of 
2D simulations as described in ISironi & Snitkovskvl j2014h . The black lines 
refer to the case of electron-positron reconnection, for comparison. The de¬ 
pendence on the numerical choice of the mass ratio mp/mc is illustrated by 
the colored curves, as indicated in the legend. 


Table 3. Results from 2D PIC simulations in electron-ion plasmas 
{niplm^ = 100) without guide field, as a function of the flow magnetiza¬ 
tion. 


O’ 

I 

3 

10 

30 

free 

0.08 

0.15 

0.22 

0.24 

(UJiU, + Ub)) 

0.52 

0.56 

0.55 

0.57 


bulk kinetic energy prior to the jet emission, and then shock waves 
are responsible for accelerating the radiating particles? 

The answer is: probably not. For shocks to be efficient in pow¬ 
ering the jet emission, we find that they need to be at least mildly 
relativistic, i.e., 7^1 > a few. To attain rough equipartition between 
t/e and Ub in the shock downstream, as inferred from modeling the 
SED of blazar sources, the magnetization in the shock upstream 
has to be in the range a ~ 0.01 - 0.1. Under these conditions, for 
most magnetic field orientations, shocks are superluminal. Charged 
particles are constrained to follow the field lines, whose orientation 
prohibits repeated crossings of the shock. So, the particles have no 
chance to undergo Fermi acceleration, which instead is required to 
explain the broadband non-thermal emission signatures of blazar 
jets. We find that the shock model can work only in a small region 
of the parameter space, with y^ei ~ 1-5 and 0i < 30°. 

If the jet remains magnetically-dominated at the dissipation 
distance (as most models of MHD jets predict), dissipation through 
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reconnection is more likely. PIC simulations of magnetic recon¬ 
nection in high-cr plasmas (see lKagan et al .120151 for a review) can 
now describe the evolution of the system for long enough times to 
probe the dynamics of the downstream region as well as the physics 
of particle acceleration. The simulations have demonstrated that 
high-cr reconnection has several attractive features to explain the 
observations. The reconnection process is fairly fast, because of 
the fragmentation of the layer into plasmoid chains. The plasmoids 
(or magnetic islands) contain the reconnection fluid and qualify as 
the emitting regions in this model. For cr > 1, the particles are 
effici ently accelerated in the co urse of the reconnection process 
(e.g.. ISironi & SnitkovskvIbOldh . In agreement with the observa¬ 
tions, the particle energy spectrum is a power-law non-thermal dis¬ 
tribution, with a harder power-law slope for higher a. In this work, 
we have shown that the reconnection process is particularly effi¬ 
cient in depositing magnetic energy into non-thermal particles, and 
at the same time that the emitting regions are in rough equipartition 
between particles and magnetic field. 

From a suite of 2D PIC simulations of relativistic reconnec¬ 
tion, we have found that, for cr > \ and zero guide field, the emit¬ 
ting regions (plasmoids) have UelUs ~ 3 and that a significant 
fraction of the the upstream energy is deposited into energetic elec¬ 
trons. More precisely, the dissipative efficiency reaches ~ 0.5 
for electron-positron plasmas (where both species can contribute to 
the emission) and pec ~ 0.2 in electron-proton plasmas. In gen¬ 
eral, the presence of a strong guide field makes the process slower 
and less efficient, and the reconnection downstream becomes more 
magnetically-dominated. Nevertheless, for a broad range of guide 
field strengths Sg/5o £ 0.5, we find that UcIUb ^ 0.3 and 
/rcc <: 0.1, which is within the range inferred from the observa- 
tions|3 This conclusions holds regardless of the flow magnetiza¬ 
tion, in the regime cr > 1 of relativistic reconnection (we tested 
from cr = 1 up to cr = 50). In fact, the measured ratio UcIUb 
might be used to directly probe the strength of the guide field in the 
reconnection region (see Figs.[3and|4j. 


The fact that cr > 1 also implies that the emitting regions may 
be characterized by relativistic motions, with outflow Lorentz fac¬ 
tors F om ~ as measured in the rest frame of the jet ( iLvubarskvl 


20051). which has nrofound imnlications for iet variabilitv (e.e.. 

Lvutikov & Blandfordll2003; 

Giannios et aLl2009l:lNalewaiko et alJ 

201 ll; Naravan & Piranll2012 

;lGiannios||2013). 


In summary, we have shown that shocks are unlikely to medi¬ 
ate the dissipation of energy in relativistic jets, mainly because the 
efficiency and subluminality constraints cannot be satisfied simulta¬ 
neously. If the bulk of the jet has cr > 1 and sufficiently small-scale 
fields prone to be dissipated, it appears that magnetic reconnec¬ 
tion can satisfy all the basic conditions for the emission: extended 
particle distributions, efficient dissipation and rough equipartition 
between particles and magnetic field in the emitting region. Dis¬ 
sipation via reconnection is a promissing process to explain the 
multi-wavelength non-thermal emission of relativistic jets. 
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